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Abstract 

In heavy-fermion compounds, the crossover from the localized to itinerant heavy-fermion state 
is observed with lowering temperature, frequently accompanied by magnetism. Ordering vectors of 
magnetism often vary with applying pressure or with substituting atoms. In Ce(Pdi_ x M x )2Al3 with 
M = Ag, Cu, the (0,0, l/2)-antiferromagnetic (AF), ferromagnetic (F), and another AF orders are 
observed for x < 0.05, 0.1 < x < 0.4, and 0.5 < x, respectively. This change in the ordering vector 
is considered to be caused by the change in the conduction-band structures. Using the anisotropic 
RKKY interaction model reflecting the spacial anisotropic distribution of the / states and also 
the conduction-band structures, we study the change in the ordering vector of Ce(Pdi_ x M :r )2Al3 
with x theoretically. As a result, the variation of the ordering vector is explained by treating the 
substitution of atoms as the conduction-electron doping, and the ordering vector of the AF state 
for x > 0.5 is considered to be (1/2, 0, 1/2). 

PACS numbers: 75.10.-b, 71.27.+a 
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I. INTRODUCTION 



Heavy-fermion systems, such as Ce-based and U-based compounds, exhibit at low tem- 
peratures the magnetic orders with various ordering vectors, which often vary with ap- 
plying pressure or with substituting atoms. In this paper, we treat the magnetism of 
Ce(Pd 1 _ a: M a; )2Al3 with M = Ag, Cu. CePd 2 Al 3 crystallizes in the PrNi 2 Al 3 -type simple 
hexagonal (SH) structure, which has the point-group and the P6/mmm space-group 
symmetry, with the lattice constants approximately a = 5.471 A and c = 4.216 A 1 . The 
magnitude of the Sommerfeld coefficient (7 ~ 380 mJmol _1 K _2 )i indicates that this com- 
pound is in the heavy-fermion state. The / states with the total angular momentum J = 5/2, 
which mainly contribute to the formation of the heavy fermion state, are divided into three 
doublets with M = ±1/2, ±3/2, and ±5/2, where M is the z (c) component of J, in the 
hexagonal crystalline electric field (CEF). CePd 2 Al 3 exhibits a large value of Xa/Xc (~ 20) 2 
in the paramagnetic phase, indicating that the CEF ground states for the / electrons are 
composed of the M = ±1/2 states mainly. In fact, the susceptibility in the paramagnetic 
phase is described well by assuming the CEF ground states as M — ±1/2, and the excited 
states as ±3/2 (33 K) and ±5/2 (800 K)2. CePd 2 Al 3 orders in the antiferromagnetic (AF) 
state of Q = (0, 0, 1/2) at Tn = 2.8 K 1 , with the saturation moment of m ~ 0.38/xb 3 - Note 
that (Q a , Qb, Q c ) denotes Q a a* + Qbb* + Q c c*, where a*, b*, and c* are the reciprocal lattice 
vectors. The change in the ordering vector is observed in Ce(Pdi_ x M x ) 2 Al 3 , in which Ag 
or Cu are substituted for Pd 4 ~ 6 . The (0, 0, l/2)-AF order region x < 0.05 is followed by the 
crossover region 0.05 < x < 0.1, the ferromagnetic (F) order region 0.1 < x < 0.4, another 
crossover region 0.4 < x < 0.5, and another AF order (the ordering vector is unreported) 
region 0.5 < x. At x = 0.2, the reversed anisotropy in the susceptibility of Xa/Xc < 1 is 
observed^, indicating that the CEF ground state has changed from the M = ±1/2 states to 
other states. Although the susceptibility cannot be fit well, it is natural to assume the CEF 
ground state for the doped system as M = ±3/2, because the estimated splitting from the 
ground state M = ±1/2 is small in the undoped system (CePd 2 Al 3 ). Moreover, the substi- 
tution of Ag or Cu gives rise to the electron doping in the conduction band. To understand 
the ordering-vector variation, we should consider the realistic conduction-band structures, 
as well as the charge distribution anisotropy of the / states. Here we report a theoretical 
study on the magnetism of Ce(Pdi_ a; M :E ) 2 Al 3 using the anisotropic RKKY interaction which 
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reflects both the spacial anisotropic distribution of the / states and the conduction-band 
structures. 

With lowering temperature (T), the / electrons in heavy-fermion compounds exhibit the 
crossover from the localized to itinerant state, as observed in the angle-resolved photoemis- 
sion spectroscopy(ARPES) for Celrln^ and UPd2Al 3 8 . For T > T* ph 50 K, the / electrons 
are localized at the energy level Ef far below the Fermi level. In the region below T*, the / 
electrons become gradually hybridized with the conduction band, to form the renormalized 
quasi-particles, namely the heavy fermions. Even in this low T region, however, the quasi- 
particles possess very large mass enhancement factor m*/m ~ z' 1 , indicating that the / 
electrons with the fraction of 1 — z ph 0.9 — 0.99 persist as the incoherent, namely localized, 
part. In this paper we discuss the relative stability of the magnetic states with different 
ordering vectors. In this case, the large fraction of the localized / electrons remaining still 
at T < T* dominantly contributes to determining the magnetism, and the rest, that is, 
the small fraction of the itinerant / electrons give only a small deviation on it. Therefore, 
the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction is responsible for the magnetism 
of these heavy-fermion compounds, even though the magnetic transition temperature, typi- 
cally 1 — 10 K, is lower than T*. The magnetic order is considered to be influenced by the 
orbital anisotropy of the / electrons, which is reflected in the susceptibility anisotropy, and 
also by the conduction-band structures. 

II. FORMALISM 

We obtain the conduction-band structures of CePd2Al3 from the band calculation using 
abinit&i^ for the non-/ counterpart system, in which La are substituted for Ce. In the band 
calculation, we adopt the Troullier-Martins pseudopotentialii, and the exchange-correlation 
term is determined according to Perdew and Wang 12 . The wave functions are expanded by 
the plane waves up to the cutoff energy of 60 Ry, and the 12 3 fc-point mesh with the special 
point technique by Monkhorst and Pack— are adopted. 

Using this conduction-band structures, we calculate the RKKY interactions between 4/ 
electrons of the Ce ions. In this case, it is necessary to take into account the effects of 
the anisotropic charge distribution of the 4/ electrons on the hybridization, as initiated by 
Coqblin and Schrieffer.— Instead of the on-site hybridizations between the conduction and 
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/ electrons (c-f hybridizations) with the plane-wave conduction states in Ref. [l4j, we adopt 
the c-f hybridizations derived from the two-center Slater-Koster integrals^ 1 ^ between the 
Ce and Pd (or Al) sites in the following way. We take account of all conduction bands 
(labeled by 7) which cross the Fermi level, where only the Al 3p, Pd 5p, and Pd Ad electrons 
are assumed to hybridize with the / electrons on the Ce sites. For the Ce 4/ electrons, we 
consider only a single ground-state CEF doublet labeled by the z component of J, M = ±|M| 
(|M| = 1/2, 3/2, or 5/2). 

Taking the number of unit cells N, the RKKY-interaction Hamiltonian is expressed as 

HrkKY = E \ M \ Ri ~ R ^ fmfiM'fjM'fjM, (1) 

ijMM' 

E m (/i, R) = —^ £b - £b , 

x a lkM e n lk i M e 

= _2_ / (g 7 fc ~ AO ~ / (fr^ ~ A*)] 

x -B 7 fcM cos(fc • R)B y k'M cos(fc' • Ji), (2) 

where /i, £ 7 fc, and /(e 7 fe — /x) are the chemical potential, the energy of the conduction band 
7, and the Fermi distribution function, respectively. B lk M consists of the terms derived 
from all process in which the 4/ electrons on a Ce sites and the conduction electron are 
scattered into the 4/ electrons on the Ce sites and conduction electrons with their angular 
momenta unchanged. The conduction electrons with given / and m are derived from various 
crystallographically inequivalent atoms, denoted as Ce, All, A12, A13, Pdl, and Pd2. B lk M 
includes two kinds of terms derived from two kinds of scattering: (1) the scattering in which 
the conduction electrons from equivalent atom are concerned; (2) the scattering in which 
the conduction electrons from inequivalent atoms are concerned. Here, we denote (1) as 
the equivalent-atom scattering and (2) as the inequivalent-atom scattering. To calculate the 
the inequivalent-atom terms in B lk M, we may be required to obtain with high accuracy the 
composition of the conduction electron using the decomposition into Wannier functions on 
the basis of the Slater-Koster method. Such a task is, however, difficult to execute. To avoid 
this difficulty, we replace the sum of such inequivalent-atom terms in B lk M by a constant Vq, 
because the sum of such terms are expected to consist of variously fc-dependent term, and 
hence to have a weak k dependence. On the other hand, for the equivalent-atom scattering, 



we calculate the corresponding terms in B^m by the use of the two-center Slater-Koster 
integrals.— 1 ^ 

Thus, we obtain the expression of B lkM as 

BjkM = 4 P~/kXnlm \ VykM — crXnlm \ 2 + V \ (3) 

Xcmlm 

where X is the label of crystallographically inequivalent atoms, a is the spin of the / electron, 
and n, I, m are the principal, azimuthal, and magnetic quantum numbers of the conduction 
electrons, respectively V^kM-aXim is the hybridization matrix element obtained from the 
Slater-Koster integrals lu-aXnim (fj)) 

aXnl m = > lM-aXnlm (Vj) t *' T * , (4) 

where j denotes all sites of X atoms, which are distant by rj from the relevant Ce site. 
The coefficient p^kXnim cannot be obtained directly from our band calculation. What we 
can obtain is only the partial density of states (DOS) of the conduction electrons derived 
from given kind of atom (Ce, Pd, or Al), with given n and I. Therefore, we assume that 
P-ykAii3 P m = -Pai3 P /9, where Pai3 P is the ratio of the partial DOS of the Al 3p electrons 
to the total conduction-electron DOS at the Fermi level, and the denominator 9 is the 
product of 3 (AU-3) and 3 (m = 1,0, —1). Similarly, we assume that p 7 fcpdi5 P m = -Ppd5p/6 
and p 7 fcpdi4dm = -Ppd4<z/10. The value obtained from our band calculation is Pai3 P = 0.27, 
Ppd5 P = 0.20, and P PdM = 0.15. 

We discuss, as an example, how to obtain |V^fcM-o-Aii3pm| • In this case, neighboring 
All atoms are distant from the relevant Ce site by t*i = a/2 + c/2, = —a/2 + c/2, 
r3 = —V\ = —a/2 — c/2, and = —V2 = a/2 — c/2. In the case that M — a — m is even, 
the Slater-Koster integrals exhibit the same value for the four distance ri.4, giving rise to 

\V lk MaA\13pm\ 2 = 16|/m-<tA11 3pm(^l) f COS 2 • fcj COS 2 (| ■ fej . (5) 

For odd M — a — m, the Slater-Koster integrals exhibit the same value for r 13 and the value 
multiplied by (—1) for r2,4, giving rise to 

|^7fcM<7A113pm| 2 = 16|/M-aA113pm(n)| 2 sin 2 ■ fcj sin 2 (| ' fc ) • (6) 

Similarly, |V^,fcAfoAi2(3)| 2 is obtained by replacing a by b (a + b) in Eqs. fl5])-(|6]), noting 
that \lM-aA\iim ( r j) | 2 is independent of i, and of j, which holds for Pd as well. Thus, 



Yl \V-ykM-aA\i3pm\ 2 is large at r (L, and H) for even (odd) M — a — m. We can give similar 

i 

discussion on Pdl (distant from Ce by r[ = 2a/3 + 6/3, r' 2 = —a/3 + 6/3, and r' 3 = 
—a/3 — 26/3) and Pd2 (distant by — r[, —r' 2 , and — r' 3 ), and obtain large \V 1 kM-aPdi5 P m\ 2 

i 

and Y2 l^7fcM-(TPdi4dm| 2 at the U and P line for M-(r-m/fl mod 3 and at the A line for 

i 

M — a — m = mod 3. 

In this way, B lkM is determined by |/Af-<7Aii3 P m| 2 , |/j\/-aPdi5 P m| 2 , and \hi-aPdi4,dm\ 2 , de- 
pendent on the ^-direction two-center Slater-Koster integrals^ 1 ^ (f per) Ah (fP 71 ) Ah (/pc)pdj 
(fP n )pd, etc. We adopt (f per) Ah (fP a )pd, (fda)pd, and V as the independent parameters, 
by assuming for simplicity that (flir)x/(flcr)x = —0.5 and (fd5)pd = 0. 

We further note that the charge distribution of the /-electron states affects the anisotropy 
of the c-f hybridizations, through the factor (7 — A Ma))/ 14 in Eq. ([3]) and the two-center 
Slater-Koster integrals. For example, the / states with M = ±1/2 have a small charge 
distribution in the basal plane, and hence a weak hybridization with the 5p and 4d states 
at Pd sites located in the basal plane, but have a strong hybridization with 3p states at 
Al sites. The resulting k dependence of B^xli shows that it has large values near the T 
point. Similarly, has large values near the T, L, and H points. On the other hand, 

the / states with M = ±5/2 strongly hybridize with the 5p and 4d states at Pd sites, and 
consequently we obtain having large values near the U and P lines. 

Since we are interested in the sufficiently low temperature region compared with the 
conduction bandwidth, / (e 7fc — /i) can be approximated by the step function 6*(/i — e 7 fc), 
where E\ M \ (/i, R) is rewritten as 

E W MR)= 4* f*- U *™*'*^'* , (7) 



/dS 
\de k /dk\ Blkm C ° S ^ ' R ^ 
f dS 

I ^^T^^lE^t^^], (8) 



Nn 6h 

"-cq 

where j £ ^ =£ dS denotes the integration over the iso-energy surface e lk = e in the first Bril- 
louin zone (FBZ). £ 7max and £ lul \ n are the maximum and minimum energies of the conduction 
band 7, respectively. fc eq are the points in the FBZ, including k itself, to which k is trans- 
lated by the symmetry operations belonging to the point group, of which the number 
of the elements Nu 6h = 24. 



Using the sine transformation <? 7 |m| (£) R) — Yl a m\M\(R) sin (mn (e — £ 7m in) where 
W 7 = e 7ma x — £ 7 min is the width of the conduction band 7, E\m\ (/^, -R) is expressed as 



E\M\ (H, -R) = ~~p2 a m|Af|(-R)«n|M|(-R) 



7 / m,n 



X 



1 A'- ^ 7 min \ j I ft £7 min 



(9) 



I m ,n(^)=Re du du . (10) 

Jo Jv U- U' 

I m ,n (y) is analytically calculated as 

sin ((m + n) z/) r . . 

Im,nt-m iy) = ; [Ci (\m\ is) - log |m| - Ci (\n\ (vr - is)) + log |n|J 



m + n 
cos ((m + n) is) 



[Si (my) + Si (n (n — 1/))] 



Im,—m 



m + n 

Si (ra/) — Si (n7r) + (— l) m+n [Si (m (7r — 1/)) — Si (rwr)] 

1 ' v. / 

m + n 

[(— l) m — 1] sin (mis) 
m 

+ vCi (\m\ u) + (tt- v) Ci (|m| (it - 1/)) - TrCi (\m\ vr) , (12) 

where Ci(x) and Si(x) are the cosine and sine integral functions. Thus, calculating nu- 
merically g-/\M\ ( £ ), which reflects the Fermi surface topology and the hybridization term 
anisotropy, we obtain the RKKY interaction. 

In the mean-field approximation (MFA), Eq. ([T]) is rewritten by the use of the Fourier 
transformations f\ M = fl M ^ ik ' Ri and f iM = ^ fkMZ~ lk ' R% as 



Qh x h 2 MM' 



— ^ K \M\Qft 1+ Q M fkiM'ft 2 M'fk2+QM, (13) 
Qfc 1 fc 2 MM' 

K W\Q = Yl E W\(^ R i ~ R j) exp[iQ • (Ri - Rj)}, (14) 

and the magnetic order parameter A\m\q for the ordering vector Q satisfies the self-consistent 
equation 

A |M|Q = -^f«ta„h%«. (15) 

Thus, the realizing ordering vector is given by Qmm, at which K\m\q takes the minimum 
K\M\Q mia , giving rise to the transition temperature T N = — K\ M \Q min / 4. Generally, if the 



RKKY interaction E\m\ (/J>, a>) between the in-plane neighboring sites is positive, the mag- 
netic order with nonzero in-plane component of Q (namely, Q a and Qb) tends to be stable. 
On the other hand, if E\ M \ (//, c) between the c-direction neighboring sites is positive, the 
order with nonzero Q c tends to be stable. 

Next, we discuss the relation between the iso-energy integrations g^\u\ (e, R) and the 
ordering vector. Since 



E\m\ (eq, R) in Eq. (J^J) is expected to increase with decreasing <7 7 |m| (eq, R) and with increas- 
ing g!yi M i (eo, R), and therefore to have a local maximum with respect to e, in the case that 
g y \M\ (e, R) changes its sign near e = e q . 

The sign of g 1 \M\ (e, R) hi Eq. (jSj) is determined mainly by the topology of the iso-energy 
surface E^k = e. For R = a, where a is the in-plane lattice vector, cos [&eq • <x] > for k 
near the A line (the central axis of the FBZ), and X] cos [^eq • a] < for fc near the outer 
edge of the FBZ. Therefore, if the iso-energy surface moves from near the A line toward 
the outer edge with increasing or decreasing e, <7 7 |m| (e, a) changes its sign, and E\m\ (e, a) 
is expected to have a local maximum. If the local maximum E\m (e m ax, a) is positive, the 
magnetic order with finite in-plane component of Q is expected to be stable near \x = £ max - 
Similarly, cos [&cq ■ c] is positive near the plane with k c = including T point (the central 

^cq 

plane in the FBZ), and is negative near the plane k c = ±1/2 including the A point (the 
top or bottom planes in the FBZ). Therefore, if the iso-energy surface moves from near 
the central plane toward the top and bottom planes with varying e, g^\M\ {e, c) change its 
sign, accompanied by a local maximum of E\m\ (e, c). If the local maximum E\m\ (^max? c) 
is positive, the magnetic order with finite Q c is expected to be stable. In the case that e is 
located near the top or bottom of the band, g 7 \M\ (e , R) is large and gL M i (e , R) is small, 
giving rise to the negative E\ M \ (e, R) which makes the F order stable. It should be noted 
that these relations between the Fermi surface topology and the ordering vector is quite 
similar to the nesting effect in the Hubbard model. 




21og2[<7 7 | M | (e )] 2 D + 



(2 -2 log 2) 
3 



(16) 
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III. RESULTS AND DISCUSSION 



We show the calculation result for (/per) ai = 0.21 eV, (fpcr)pd = 0.21 eV, (fda)pd = 0.14 
eV, Vq = 0.17 eV, and the depth of the bare / level Ef = —2.5 eV. Among the conduction 
bands which we obtained, the 16th and 17th bands possess the Fermi surfaces. The density 
of states (DOS) of the 15-18th bands, which are considered in the calculation of the RKKY 
interaction, are plotted in Fig. [TJ It follows that CePd2Al 3 (the band structures with jj, — 0) 
exists in the crossover from the 16th-band dominant region to the 17th-band dominant one. 
Several of the iso-energy surfaces for the 16th and 17th bands are plotted in Figs. Hfa)-(d) 
and E^a)-(d), respectively. From e = —0.6 eV to 0.0 eV, the iso-energy surface of the 16th 
band approaches the central a*6*-plane (k = (k a , kb, 0)), giving rise to the change in the sign 
of #16 1/2 (e, c) from negative to positive at e ~ —0.5 eV as plotted in Fig. 0] (a). On the other 
hand, the iso-energy surface of the 17th band moves from the vicinity of the A line (the 
central axis with k = (0,0, k c )) to outside from e = 0.0 eV to 0.6 eV, and simultaneously 
approaches the central a*6*-plane (k = (k a ,kb,0)) at 0.4 eV < e, giving rise to the change 
in the sign of #171/2 (e, a) from positive to negative at e ~ 0.4 eV and that of #171/2 (e, c) 
from negative to positive at e ~ 0.5 eV, as plotted in Figs. @] (b). Because (e, R) is not 
largely different from <7 7 i/2 (e,R), it is not plotted here. 

These behaviors of g^\M\ (e> R) give rise to the local maximum of E\m\ (e> c ) near —0.3 
eV and near 0.6 eV, and that of E\m\ (e, a) near 0.3 eV, as plotted in Figs. |5] (a) and 
(b), respectively. Thus, the (0,0, l/2)-AF and (1/2, 0, l/2)-AF states become stable near 
[i = —0.2 eV and near 0.3 eV, as plotted in Figs. E] (a) and (b). In the intermediate region 
H ~ 0, where \i is located at the crossover from the 16th dominant region to 17th band 
dominant one, the F order is stable. 

Next, we discuss the dependence of the phase diagram obtained here on the conduction- 
band structures and its consistency with the observed ordering-vector variation in 
Ce(Pdi_ I M a; )2Al3. The obtained band structures with \i = correspond to the conduction- 
band structures of CePd2Al 3 . It is difficult to calculate the conduction-band structures of 
the mixed crystal Ce(Pdi_ x M :r )2Al3, and hence we describe simply the conduction band of 
Ce(Pdi_ x M2;)2Al3 as the system to which 2x electrons per unit cell are doped, by which 
increases. The calculated phase diagrams for M = ±1/2 and ±3/2 are plotted versus \x 
in Figs. El (a) and (b), respectively. For M = ±1/2, the (0, 0, l/2)-AF, F, (1/2, 0, l/2)-AF, 
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16th 




FIG. 1: (Color online) The density of states (DOS) versus e. The thin solid, thick solid, thick 
dashed, and thin dashed lines denote the DOS of the 15th, 16th, 17th, and 18th bands, respectively, 
e = corresponds to the Fermi level for CePd2Al3. 




FIG. 2: (Color online) The iso-energy surfaces of the 16th band for (a): e = —0.6 eV; (b): e = —0.4 
eV; (c): —0.2 eV; (d): 0.0 eV (calculated Fermi surface). The T point is located at the center. 

(1/3, 1/3, 0)-AF, and (0,0,l/2)-AF orders are most stable in the regions —0.31 eV < /i < 
-0.09 eV, -0.08 eV < /i < 0.13 eV, 0.20 eV < fi < 0.35 eV, 0.42 eV < /i < 0.51 eV, and 
0.56 eV < n < 1.03 eV, respectively. For M = ±3/2, the (0, 0, l/2)-AF, F, (1/2, 0, l/2)-AF, 
(1/2, 0, 0)-AF, and (0,0, l/2)-AF orders are most stable in the regions —0.31 eV < /i < 
-0.07 eV, -0.06 eV < fi < 0.13 eV, 0.17 eV < fi < 0.41 eV, 0.42 eV < /i < 0.53 eV, 
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FIG. 3: (Color online) The iso-energy surfaces of the 17th band for (a): e = 0.0 eV (calculated 
Fermi surface); (b): 0.2 eV; (c): 0.4 eV; (d): 0.6 eV. 
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FIG. 4: (Color online) (a): g W i/2 (e, R) versus e; (b): g 171 / 2 (e ) -R) versus e. The solid and the 
dashed lines denote the R = a and R = c, respectively. 

and 0.56 eV < fi < 0.99 eV, respectively. The difference in the phase diagram between 
M = ±1/2 and ±3/2 are not significant, except for the region fi ~ 0.5 eV, suggesting that 
the difference in the charge distribution anisotropy between M = ±1/2 and 3/2 states is not 
crucial in the determination of the ordering vector, compared with the shift of fi caused by 
the substitution of atoms. We note further that the sequence of the stable phase is nearly 
invariant, except for the region /i ~ 0.5 eV, even if is assumed as constant, indicating 
that the k dependence of is not important unless it is too strong. 

We plot the range of the stable ordering vectors versus x for M = ±1/2 and ±3/2, in 
comparison with the experimental results^-, in Fig. [71 We obtain the sequence of the AF 
order with Q = (0, 0, 1/2), the F order with Q = (0, 0, 0), and another AF order, of which 
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FIG. 5: (Color online) (a): EWi (e,a) versus /i; (b): E\ M \ (e,c) versus fi. The solid and the dashed 
lines denote E\ M \ for the M = ±1/2 and ±3/2 states, respectively. 
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FIG. 6: (Color online) /x-T phase diagram for (a): M = ±1/2; (b): M = ±3/2. 
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FIG. 7: (Color online) The range of the experimentally^— (top) and theoretically (M = ±1/2 in 
the middle, and ±3/2 in the bottom) obtained stable ordering vectors versus x. 
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the ordering vector is (1/2,0, 1/2). This sequence is in good agreement with the observed 
one for Ce(Pdi_ x M.j;)2Al3 with increasing x. In order for the theoretical results to correspond 
with the experimental ones quantitatively, it is required to shift x by approximately +0.1, 
which corresponds to the change of the conduction-electron number by 0.2. 

To dissolve this discrepancy, it must be allowed to assume a small change in the 
conduction-band structures. If the 18th conduction band, which are thoroughly unoccu- 
pied in our band calculation as plotted in Fig. [TJ are lowered so that it possesses small 
electron Fermi surfaces with the filled electron number ~ 0.2, we may obtain effectively 
smaller electron filling numbers for the 16th and 17th bands. Under this modification in 
the conduction-band structures, we can relate the (0, 0, l/2)-AF order for fi ~ —0.1 eV, the 
F order in the vicinity of \i — 0, and the (1/2,0, l/2)-AF order for fi ~ 0.2 eV to the AF 
order of CePd 2 Al 3 , the F order of Ce^di-zM^Ahj for 0.1 < x < 0.4, and the AF order for 
x > 0.5, respectively. The (0,0, l/2)-AF order is caused by the RKKY interaction through 
the structure of the 16th conduction band, while the (1/2, 0, l/2)-AF order is due to the 
17th band, and the F order is caused by the crossover from the 16th to the 17th band. 

Finally, we discuss the justification of the mean-field approximation (MFA). In our model, 
the stable magnetic order is determined mainly by the signs of the RKKY interaction be- 
tween the neighboring Ce sites. Therefore, the stable ordering vector will not alter substan- 
tially even if the approximation is improved beyond the MFA, but the transition temperature 
will be largely changed. As a conclusion, the MFA is sufficient inasmuch as the relative sta- 
bility between the various magnetism is concerned. A calculation based on an improved 
approximation is left in a future issue. 

IV. CONCLUSION 

In this paper, we have studied the variation of the ordering vectors of Ce^di-zM^Ala on 
the basis of the anisotropic RKKY interaction model, using the realistic conduction-band 
structures obtained by the band calculation. It is understood that the conduction-band 
structures play an important role in determining the ordering vectors, in the way similar to 
the nesting effect in the Hubbard model. By treating the substitution of Ag or Cu for Pd 
atoms as the electron doping to the conduction bands of CePd2Al 3 , which are obtained from 
the band calculation, we have obtained the variation of the ordering vectors consistent with 
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experimental results for Ce(Pdi_ x M x .) 2 Al 3 . The (0,0, l/2)-AF order of CePd2Al 3 is caused 
by the movement of the iso-energy surface of the 16th band to the central plane of the FBZ 
with increasing e. On the other hand, the AF order for x > 0.5 is caused by the movement 
of the iso-energy surface of the 17th band, for which the ordering vector is predicted to 
be (1/2,0,1/2). The F order for 0.1 < x < 0.4 is caused by the crossover of the main 
conduction band from the 16th to the 17th band. Thus, we have obtained the variation of 
the ordering vectors with increasing /x, which is consistent with that for Ce(Pdi_ x M x )2Al3, 
using the anisotropic RKKY model with the realistic conduction-band structures. 
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